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Within QRPA we achieve partial restoration of the isospin symmetry and hence fulfillment of the 
requirement that the 2i//3/3 Fermi matrix element MJ." vanishes, as it should, unlike in the previous 
version of the method. This is accomplished by separating the renormalization parameter g^j, of 
the particle-particle proton-neutron interaction into the isovector and isoscalar parts. The isovector 
parameter need to be chosen to be essentially equal to the pairing constant gpair, so no new 

parameter is needed. For the Ou/3/3 decay the Fermi matrix element M^" is substantially reduced, 
while the full matrix element M^'^ is reduced by ~ 10%. We argue that this more consistent approach 
should be used from now on in the proton-neutron QRPA and in analogous methods. 
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I. INTRODUCTION 

Answering the questions whether total lepton number is a 
conserved quantity or not, and thus whether neutrinos are 
massive Majorana fermions, is a crucial part of the search for 
the "Physics Beyond the Standard Model". Consequently, 
experimental searches for the Oi//3/3 decay are pursued world- 
wide (for a recent review of the field, see e.g. [T]). However, 
interpreting existing results and planning new experiments 
is impossible without the knowledge of the corresponding 
nuclear matrix elements. 

The nuclear matrix elements A/"" of the 0i^/3(3 decay must 
be determined using nuclear structure theory, and the choice 
of the appropriate approximations is a crucial part of that 
task. Some of the methods employed for evaluation of the 
M'^'^ , in particular those that begin with the transforma- 
tion from particles to quasiparticles to account for the like- 
nucleon pairing ( see e.g. [MS])) use wave functions that do 
not exactly conserve the particle number. The number of 
protons and neutrons is usually conserved on average or, in 
some cases, it is restored by the particle number projection. 
In either case, until now no attempt was made to check that 
the isospin, which is known to be, to a very good approxi- 
mation, valid quantum number in nuclei, remains as such in 
the resulting wave functions that are obtained by solving the 
corresponding equations of motion. 

It is well known that by the proper treatment of the quasi- 



particle interaction the broken symmetries can be restored. 
Naturally, exact calculation would restore the broken sym- 
metries exactly. However, even with the approximate, RPA- 
like treatment, it is possible to partially restore some of the 
broken symmetries. In this work we show, following basi- 
cally the suggestions made initially in Ref.[n], how this can 
be done in the case of isospin and by doing that the values 
of the Fermi nuclear matrix elements, both for the 2i^/3/3 and 
0iy(3/3 decays, are substantially modified. Even though the re- 
sulting total M"'^ nuclear matrix elements are changed only 
by 10%, it is worthwhile, and certainly more consistent, 
to use in future the prescriptions described below. 



II. FORMALISM 

Assuming that the Oz//3/3 decay is caused by the exchange 
of light Majorana neutrinos, the half-life and the nuclear ma- 
trix element are related through 
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^1/2 



(1) 



where G'^'^{Q, Z) is the calculable phase space factor, {mpis) 
is the effective neutrino Majorana mass whose determination 
is the ultimate goal of the experiments, and M'^'^ is the nu- 
clear matrix element consisting of Gamow- Teller, Fermi and 
Tensor parts. 
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the 



MGT(l-Xi=^M + XT) , (2) 



matrix element 



ratios xf 



where xf and xt 

M^^/M^^j, and xt = M^'^/M^^. (In the literature a dif- 
ferent notation is sometimes used, xf — Hs^^gt ) 
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The main GT part, M^"^, can be somewhat symbohcaUy 
written as 



M'gt 



(/ISjfeO-/ • o-kT^^T^ H{rik,E)\i) , 



(3) 



where H(rik, E) is the neutrino potential described in detail 
in [S] and r/^ is the relative distance between the two neu- 
trons that are transformed in the decay into the two protons. 
Analogously, the Fermi matrix element is 



U\-i:ikT+T+H{rik,E)\i) 



(4) 



Note that these Qvfiji matrix elements are expressed in the 
closure approximation; its applicability is also discussed in 
[5]. However, the results reported later in this work were ob- 
tained without using the closure; instead explicit summation 
over all virtual intermediate states was performed. 

The half- life of the experimentally well studied 2i//3/3 decay 
depends formally on two nuclear matrix elements 
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G2'^(g,Z)[Af, 
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The Gamow- Teller 2v(ij3 matrix element is 



(/||i;feCTfeT+||m)(m||i;;0-;T+||i) 



En 



{Kh + Mf)/2 



(5) 



(6) 



where the summation extends over all 1+ virtual intermedi- 
ate states. In that case the closure approximation is not a 
valid approach but can be formally introduced by defining 
the corresponding closure matrix element M'^l' when replac- 
ing the energies Em by the proper average value E2u- Thus, 



f2v 



M'^ X {E^^^GT - (M, + Mf)/2) 



(7) 



Formally, in the description of the 2i>/3/3 decay also appears 
the Fermi matrix element 



Mir"" 



(/||S]fcr+||m)(m||I],r+||z) 
' £;,„-(Af, + M/)/2 



(8) 



where the summation extends over all 0"'' virtual intermedi- 
ate states, and its closure form is 
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Fid) 



(/|S/.r+T+| 



M 



F(ci) = X (£;2.-i. - (M, -I- Mf)/2) 



(9) 



The ground state \i) of the initial nucleus has isospin T = 
T:, = {N - Z) /2 while the final state |/) has isospin T-2 = 
Tz — {N — Z — 4)/2. Since the operator SfcT^f just changes 
the isospin projection and cannot change the total isospin, 
it is obvious that when isospin is a good quantum number 
both Fermi matrix elements must vanish. 
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M 



Fid) 



0, 



(10) 



since the average energy denominators in Eq. ^ are nonva- 
nishing. 



Until now, within QRPA, PHFB, EDF and IBM-2 meth- 
ods the validity of condition Eq.(10l has not been usually 



tested (ISIH])- When resuhs were published, Mf and M|.'^^;) 
do not vanish and are, in fact, comparable to and 
^'^GT{ci)' respectively. Despite that, when evaluating the 
2i//3/3 half-life the Fermi matrix element was usually simply 
neglected. 

As we show further, in the usual application of QRPA the 
condition Eq.(lO) is not obeyed. Instead, the magnitude of 
MJ." is numerically comparable to the magnitude of Mq^ 
as just pointed out. In addition, for the Qv^jS decay, within 
QRPA the ratio xf ~ —0.5 while in the nuclear shell model, 
where isospin is a good quantum number by construction, the 
condition Eq.(10l is, naturally, obeyed and xf ~ —(0.2 — 0.3) 

Where does this problem in QRPA method originate? The 
method begins with the Bogoliubov transformation relat- 
ing the particle creation and annihilation operators aj^, cijm 
with the quasiparticle creation and annihilation operators 
Cjjy^jCjm- By solving the BCS equations one includes the 
neutron-neutron and proton-proton isovector pairing inter- 
actions. 

At this stage several symmetries are broken. The numbers 
of protons Z and neutrons are no longer exact, but valid 
only on average. In addition, since the neutron-proton part 
of the isovector pairing interaction is neglected, additional 
source of isospin violation is introduced. It turns out that it 
is relatively easy to remedy this additional effect and restore 
the isospin conservation, at least in part, as explained further 
here. As the RPA ( and QRPA) is derived from the equation 
of motion for bifermionic operators (treated in the quasibo- 
son approximation), symmetries of the model Hamiltonian 
can naturally be fulfilled in that approximation. 

To proceed further, the equations of motion need to 
be solved. Within the QRPA method the forward- and 
backward-going amplitudes X and Y that are needed for the 
evaluation of the nuclear matrix elements, as well as the cor- 
responding energy eigenvalues uj^ , are determined by solving 
the eigenvalue equations of motion for each angular momen- 
tum and parity J'^ 



A B 

-B -A 



The matrices A and B are (see e.g. [11] ) 



pn,p'n' 

(0|[(c,c„)(^^^)\[7?,(ct,ci,)^'^*'^]]|0) 

^pp'^nn' {Ep -\- -E^ii) 
{UpVnUp'Vn' + VpU„Vp^U„') X 

2gph{pn~\j\V\p'n'-\j} 

) X 

'^9pp{pn,J\V\p'n',J) , 



(11) 



(12) 



and 



B 



pn,p'n' 



(13) 



3 



= ~{UpV„Vp'Un' + VpUnUp'Vn') 

x2gph{pn~\j\V\p'n'-\j) 

x2gpp{pn,J\V\p'n',J) , 

where Ep, are the quasiparticle energies. 

The definitions, Eqs. (12) and (13), contain two renormal- 
ization adjustable parameters gph for the particle-hole inter- 
action, and gpp for the particle-particle interaction. While 
9ph = 1-0 is typically used, it is customary to adjust gpp so 
that the experimentally known half-life of the 2i//3/3 decay is 
correctly reproduced [2]. But the particle-particle neutron- 
proton interaction governed by gpp actually contains two 
kinds of interaction matrix elements, isovector and isoscalar. 
Thus, to be consistent with the treatment of the like particle 
pairing, one should separate the T = 1 part from the T = 
part, i.e. replace 



gpp{pn,J\V\p'n',J) 
gl=\pn,J,T=l\V\p'n\J,T 



1) 



(14) 



+<7j=°(pn, J, T = Q\V\p'n', J,T = 0) 



and adjust the parameters gj^^^ and gp^'^ independently. 
To part ially restore the isospin symmetry and achieve that 
Eq.(10l is obeyed, it is sufficient to choose gpj^^ ~ gpair- 
(Tharthe coupling constant of the isovector proton-neutron 
particle-particle force should be close, or identical, to the 
pairing strength constant, was recognized already in the early 
works on the QRPA application to the /3/? decay that used a 
schematic, (5-force interaction, see Ref. [T2]). 



III. DETERMINATION OF THE PARAMETER g'^=^ 

When solving the BCS pairing equations, it is customary 
to slightly renormalizc the strength of the pairing part of 
the realistic nucleon-nucleon interaction so that experimental 
pairing gaps are correctly reproduced. Thus, four adjusted 
parameters {dpp,dl/J ,dn}i,dn2) are introduced (see, e.g. [2- 
[S]) representing the adjustments needed to describe the neu- 
tron and proton pairing gaps in the initial and final nuclei. 
The values of these parameters as well as their averages for 
selected /3/3-decay candidate nuclei are displayed in Table |Tj 
(The Table entries are for two variants of the nucleon-nucleon 
interaction and one choice, of large size (21/23 levels, oscil- 



lator shells iV = — 5 with the addition of the 



■13/2 



and 



■11/2 



for the nuclei heavier than ^^^Sn), of the single particle level 
scheme. The results for other choices are not very different.) 
In several cases in Table |l] we encounter magic numbers of 
neutrons or protons. In those cases the BCS treatment is 
inappropriate and hence the corresponding entries are miss- 
ing there. For the case of ^*Ca we considered two variants. 
In the listed one we assumed that there is no pairing in the 
doubly magic ^®Ca. In the other variant we assumed that 
the values Ap = 2.18 MeV and A„ = 1.68 obtained from 



the usual odd-even mass difference with the five point for- 
mula represent the pairing gaps; the resulting gpj^^ is rather 
similar to the values listed in Table HI 
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FIG. 1: (Color online) Dependence ol the 21/13(3 matrix elements 
MJ-" and Afgy on the isovector coupling constant Qp^^ ■ This 
example is for ^^Ge and ^'^"Te, the Argonne V18 potential, and 
the isoscalar coupling constants gjp"" = 0.750 and 0.783. 



The example in Fig. [Tjshows how the matrix elements Mp'^ 
and Mq^ behave when the isovector coupling constant g^^^ 



r=i 
pp 



is changed while the isoscalar gp~ is kept constant. As one 
can see, the Fermi matrix element Mp^ decreases and crosses 
zero, with increasing g^^^ , while the Gamow- Teller matrix 
element remains constant. This is a typical case, and we 
can now choose g'^^^ such that M'^_^i and hence also M'p' 
vanish. Those values of qE^^ are shown in the last column 
of Table H 

It follows from the entries in Table [J that the renormaliza- 
tion parameter g^^^ that assures the validity of the Eq. ( 10 ) 
is indeed very close to the average d = gpair of the pairing pa- 
rameters d\/ . In few rare cases, in particular in semi-magic 
nuclei, the difference between gp^^ and d is ~10% (but not 
more). As shown in Ref. |9] the ratio gpj^^ /d remains essen- 
tially unchanged when the size of the single-particle basis is 
modified. 

In this work, our choice is to renormalize gpj^^ indepen- 
dently, but very close, to d. On the other hand, it turns 
out that the M^'^ depends sensitively only on g'^^^ , so that 
this parameter can be still adjusted such that the half-life 
of the 2J//3/3 decay is correctly reproduced, exactly as done 
before. In fact, the previously used common value of gpp, 
and the new parameter g^^'^ , are essentially the same as we 
demonstrate in the next Section. 
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TABLE I: Renormalization parameters of the pairing interaction, their average and the T — 1 renormalization constant gjp ^ adjusted 
such that Mpl^^i-^ and MJ." vanish. 



nucleus 


NN pot. 


number 






Ujpp 


) 

^nn 


d 


9pp 






of s. 


p. lev 














^*Ca 


Argonne 


21 


lev. 


— 





1.075 


0.988 


1.034 


1.031 




CD-Bonn 


21 


lev. 


— 





0.985 


0.903 


0.944 


0.944 




Argonne 


21 


lev. 


0.930 


1.074 


0.970 


1.106 


1.020 


1.038 




CD-Bonn 


21 


lev. 


0.863 


0.983 


0.899 


1.013 


0.940 


0.958 




Argonne 


21 


lev. 


0.869 


1.085 


0.930 


1.131 


1.004 


1.032 




CD-Bonn 


21 


lev. 


0.808 


0.995 


0.864 


1.038 


0.926 


0.955 


^^Zr 


Argonne 


21 


lev. 


0.923 


0.768 


1.000 


0.962 


0.913 


0.984 




CD-Bonn 


21 


lev. 


0.856 


0.704 


0.926 


0.881 


0.842 


0.907 


^°°Mo 


Argonne 


21 


lev. 


1.019 


0.960 


1.041 


0.979 


1.000 


1.008 




CD-Bonn 


21 


lev. 


0.946 


0.883 


0.966 


0.900 


0.924 


0.933 




Argonne 


21 


lev. 


1.000 


0.975 


1.025 


0.945 


0.986 


0.979 




CD-Bonn 


21 


lev. 


0.930 


0.895 


0.954 


0.871 


0.913 


0.908 




Argonne 


21 


lev. 


1.017 


0.971 


— 


0.919 


0.969 


0.922 




CD-Bonn 


21 


lev. 


0.949 


0.895 


— 


0.847 


0.897 


0.852 




Argonne 


23 


lev. 


— 


1.001 


0.929 


1.000 


0.977 


0.988 




CD-Bonn 


23 


lev. 


— 


0.918 


0.860 


0.917 


0.898 


0.913 


128yg 


Argonne 


23 


lev. 


0.881 


0.968 


0.926 


0.999 


0.944 


0.988 




CD-Bonn 


23 


lev. 


0.816 


0.889 


0.857 


0.918 


0.870 


0.914 


130yg 


Argonne 


23 


lev. 


0.845 


0.970 


0.920 


1.000 


0.934 


0.989 




CD-Bonn 


23 


lev. 


0.783 


0.891 


0.852 


0.918 


0.861 


0.915 




Argonne 


23 


lev. 


0.851 


0.912 


0.917 


0.963 


0.911 


0.973 




CD-Bonn 


23 


lev. 


0.790 


0.840 


0.850 


0.887 


0.842 


0.903 




Argonne 


23 


lev. 


0.782 




0.885 


0.926 


0.864 


0.950 




CD-Bonn 


23 


lev. 


0.726 




0.821 


0.853 


0.800 


0.881 



IV. RESULTS AND DISCUSSION 

In the previous Section we explained how the parame- 
ter gpj^^ is determined. The determination of the other 
renormahzation parameter gj^^" is analogous and follows 
the suggestion made long time ago in Ref. [5] as already 
stated. We fit 5^""° from the requirement that the calcu- 
lated values of the full 2i//3/3 matrix elements Mq^ agrees 
with their experimental values. For most nuclei in Table I 
the half-lives T^j2 have been measured; we use the recom- 
mended values in Ref. [14] plus the ^'^^Xe half-life of Refs. 
[T51 [TB] . But for several nuclei in that Table, and in the 
Tables that follow, the half-life remains unknown. In those 
cases we proceed as follows: for ^^°Pd we use the estimate 
of Ref. [Ilj that uses the single-state dominance assump- 
tion, for ^^''Sn and ^'^''Xe we use an interval of possible M^^j, 
values < M^^j, < 0.2(0.1)MeV-i for i24Sn(i34Xe) , respec- 
tively. For these two nuclei we show in Tables II - IV the 
results with gA = 1.0 for the upper limit of Mq^, and with 
qa = 1-27 for the lower limit M^:^ — 0. Our results for these 
two nuclei therefore reflect our incomplete knowledge of the 



corresponding 2i> half-life. 

Before presenting the results for the Oz//3/3 nuclear matrix 
elements, several comments are in order. Since the main 
effect considered here is the change in Mp'^ and the associated 
change in Mp^, lets analyze these changes using the radial 
dependence of M'^_^i explained in Ref. [5]. In Fig. [2] the 
functions C'^_^i(r) with the old and new parametrization of 
gpp are plotted, together with the function C^rp_^i{r) (scaled 
by 1/3 for clarity). As one can see, with the new ffjp^^ the tail 
of C|?l^;(r) becomes more negative and therefore its integral 
vanishes, as required. Let us remind ourselves that 

/■OO 

M'/^ci = / Cf_,i{r)dr , (15) 
Jo 

and in analogy 

Ml^T-ci = / Cl^T-ci{r)dT . (16) 

Another comment concerns the fact that, as we will see, 
with the new parametrization xf ~ —(0.3 — 0.4), or more 



5 



0.6 



I=F (old par.) 

I=GT (old/new par.) 

I=F (new par.) 




relations must be modified: 



FIG. 2: (Color online) Functions Cp'^{r) with old and new 
parametrization and, for comparison, the function Cgy(r) (scaled 
by 1/3) for clarity) is also shown. This is the case of ^^Ge. 



precisely xf ~ —1/3. Somewhat similar conclusion is ob- 
tained in the shell model [TO] where the isospin is conserved 
by definition. The explanation is based on the fact that the 
ground states of even-even nuclei consist dominantly of the 
= 0+,r = 1 Cooper pairs that, in turn, are mostly in 
the S — 0, L — state. Since such states are eigenstates of 
the operator ai ■ a-i with eigenvalue -3, our conclusion simply 
follows. 

In Fig. [3] we show examples of the decomposition of the 
function C^1^;(r) into their 5 = and S = \ components. 
These are rather typical cases. The dominance of the 5 = 
component in the pure pairing case (the upper panel in Fig. 
|3]) is easily understood. However, that feature is still present 
in the realistic case with g^-p ^ 0.0, hence our finding that, 
usually, Xf ~ —1/3. 

For both modes, Qvfip and 2i/(3(], we can find relations 
between the Fermi and Gamow- Teller parts and their 5 = 
and 5=1 components. These relations are exact in the 
closure approximation and when the higher order weak cur- 
rents ( and thus the tensor part Mjl^ ) are neglected and 
the nucleon form factors have the same cut-off values for the 
vector and axial vector parts. In addition, since the neu- 
trino potentials indirectly depend on the assumed averaged 
energy, these E values must be chosen to be the same for 
the Fermi and Gamow- Teller matrix elements. Using the 
properties of the cti • tT2 operator and that Mj," = with 
our new parametrization we find that only one of the four 
components is independent and 



M^r-pis = 0) 



-M^r-Fis^i) 



(17) 



- -M'^i-GTiS = 1) = -M^r^ariS = 0)/3 . 
For the Ov mode, however, Mp" ^ and hence the above 



M^''(5 = 0) 



A/0'-'(5= 1) : 



0). 



0)/3 



(18) 



Two components are independent in this case. In realistic 
case these relations are not exact, but still valid in a reason- 
able approximation. 

We will return to the discussion of the xf values obtained 
by different approximate methods in the next Section. 



1.2 R 
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FIG. 3: (Color online) In the upper panel (a) the function 
C^_^i{r) is shown for the pure pairing case, i.e. g^^^ — QpjT^ = 
dph ~ 0.0, separated into the S — and S = 1 components. In 
the lower panel (b) the function C^_^i{r) is shown for (?^^ = 
1.038 and ^j^" = 0.750 again separated into its 5* = and S = 1 
parts. The sum function is also displayed. The dominance of the 
5 = component is clearly visible in the upper panel. In the 
lower panel the two components when integrated over r are, nat- 
urally, equal and opposite. The 5 = part, however, clearly is 
considerably larger in absolute value than the 5=1 part, at all 
r values. This is the case of ^^Ge. 
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In tables |II] and |III| we compare the resulting matrix ele- 

th it 

using the old parametrization {g^^^ 



ments Mj,", M'^j. and M^'^ with its components evaluated 



dpp) with the 
"^=0 is fitted 



T=0 

^pp 9pp 

new results, where g^^^ « gpair and where 
to the known experimental values of Mq^. The calculations 
were performed for the unquenched value = 1.27 as well 
as for gA ^ 1.0. The quantities M'^"" = M°'' x (g^/1.27)2 
are also shown, as well as xf- Calculations in both tables 
were performed within the standard QRPA with all the usual 
ingredients, i.e. including the higher order weak currents, nu- 
cleon form factors, and the short-range correlation treatment 
ofRef. pj]. 

Lets explain briefly again (for more details see Ref. [5]) 
how the quenching is taken into account with our method. 
Since we adjust the isoscalar particle-particle renormaliza- 
tion constant gpj^'~' in such a way that the experimental half- 
life of the 21/(30 is correctly reproduced, by changing the 
effective value of the axial current coupling constant g^ we 
are forced to change also the parameter gpj^^, albeit only 
slightly. Those changes are visible in the third columns of 
Tables II and III. Since with smaller gA the parameter gpj^*^ 
slightly decreases, the corresponding Mq'j, matrix element 
increases. However, the Ov/3f3 decay rate, proportional to 
the (M'^")^, naturally, decreases. 

In that context it is worthwhile to point out another fea- 
ture of the new parametrization. The Fermi matrix element 
Mp" is associated with the weak vector current, and as such 
should not be affected by the axial current quenching. With 
the old parametrization, with a single gpp, that was not 
quite true, as seen in Tables II and III. However, with the 
new parametrization where isospin symmetry is partially re- 
stored, the Mp" becomes independent of the effective value 
of gA as it should be. (The tiny changes in Tables II and III 
are round-off errors.) 

From the tables one can see that the new parametrization, 
leading to Mp" — 0.0, leads to a substantial reduction of 
the M^" component of M^" and overall 10-20% reduction 
of the final Af^" nuclear matrix elements. It is encouraging 
that both variants of the M^"^ matrix elements for ''^Ca are 
now rather close to the results of nuclear shell model eval- 
uation. (With gA — 1.27 our M'^" values are 0.54 in the 
listed case and 0.71 in the variant where the even-odd mass 
differences are treated as arising from pairing, both with the 
Argonne V18 potential and 0.59 (0.77) with the CD-Bonn 
potential , while the shell model values are 0.59 in Ref. [TH] 
and 0.82(0.90) for the Argonne VIS (CD Bonn) potential in 
Ref. [El .) Note that only in the case of '*®Ca the full oscillator 
pf shell is included and hence the Ikeda sum rule is fulfilled 
in the nuclear shell model treatment. We are, naturally, well 
aware of the fact that to apply QRPA in the case of ''^Ca 
is questionable; our results should be treated with that in 
mind. 

Finally, in order to better visualize effect of the new 
parametrization of the particle-particle interaction, we show 
in Fig. [4] an example of the multipole decomposition of the 
matrix element Mpy. One can see there that the contribu- 



tion of the intermediate multipole 0^ is drastically reduced 
with our choice of gpj^^, while all the other multipoles are 
affected only slightly or not at all. This is, in some sense, 
analogous to the situation with Mq^ where the parameter 
gjp^^ affects mostly the intermediate 1"*" states, while all the 
other multipolarities are affected much less. 



I oldparam. 
I new param. 




FIG. 4: (Color online) Multipole decomposition of the matrix 
element M^" . The results with the old and new parametrizations 
are compared. Note the dominant effect for the 0"'' multipole, and 
the relatively small effects for the other multipoles. This is the 
case of ^^Ge. 

We compare in Fig. [5]the M*^" matrix elements for all con- 
sidered nuclei evaluated with the old and new parametriza- 
tions of gpp. The smaller values of M^" in "SCa, ^^^Cd, 
^^"^Sn, ^^^Xe and to some extent also in ^^Zr are related to 
the magic or scmimagic nucleon number in these nuclei, and 
thus to the reduced pairing correlations in them. 



V. COMPARISON OF THE xf VALUES 
EVALUATED BY DIFFERENT METHODS 

As we argued in this work, the result of the new 
parametrization of the particle-particle interaction, that par- 
tially restores the isospin symmetry and leads to the correct 
Mp" = value, is the reduction of the Fermi part Mp'^ of the 
Oiy/SfS nuclear matrix element. At the same time, the largest 
component of that matrix element, Mq^, remains essentially 
unaffected. One can see that most clearly by considering the 
quantity Xi=^, the ratio M^'^/Mg^. 

In Table [TV] we compare the xf values obtained with dif- 
ferent methods. (Analogous table, naturally without our new 
results, appears in Ref. ^0^ in their Table VII. However, as 
we already mentioned, their definition of xf c ont ains an ex- 
tra factor (gv/gA)'^-) One can see in the Table IV that in the 
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FIG. 5: (Color online) Nuclear matrix elements M''" evalu- 
ated with the new parametrization developed in this work (filled 
squares) compared with the old method {gpjT^ — Qpp^ = Qpp) 
(empty circles). This is QRPA with qa = 1.27 and large size 
single particle level scheme, as in Table I, evaluation using the 
Argonne V18 potential. 

nuclear shell model, and in our QRPA calculation with the 
new parametrization of gpp, the xf values are substantially 
smaller than in the previous approaches. (In IBM-2 the xf 
are very small when neutrons and protons are in different 
shells. That is an artifact of the model where only one shell 
in each system is included.) 

In the shell model, and in our new QRPA calculations, the 
Xf values are relatively close to -1/3, the value one would 
obtain in pure 5 = states. However, in the shell model 
the Xf values are systematically smaller than in our version 
of QRPA. Why this is so remains to be understood. (To 
be really precise, xf = ^1/3 would arise for pure 5 = 
when the higher order terms in the weak current are absent, 
when in the nucleon form factor the cut-off parameters for 
the vector and axial vector currents are the same and the 
average energies E are the chosen to be the same in both 
neutrino potentials.) As we pointed out before, while the 
5 = component is large, the other parts, in particular 
5=1, are clearly present. 

We may notice that the QRPA values of xf are always 



smaller with the quenched value qa — 1.0 compared to the 
unquenched value gA = 1.27. That trend continues when 
the amount of quenching is increased, e.g. to gA = 0.8 
where xf values are really quite close to -1/3. However, the 
question of quenching of the QvjSp matrix elements remains 
open, and in particular how to treat it properly in the 
QRPA goes beyond the scope of the present paper. 

VI. CONCLUSIONS 

By separating the particle-particle neutron-proton inter- 
action into its isovector and isoscalar parts, and renormal- 
izing them each separately with its own fitted parameters 
9pp^^ and g^^^ , we have achieved the partial restoration of 
the isospin symmetry and fulfillment of the requirement that 
M'P — 0.0. This has been done essentially without introduc- 
ing new parameters, since g^^"^ ~ gpair as required by the 
isospin symmetry of the particle-particle force. At the same 
time the isoscalar parameter g^^^ is fitted from the require- 
ment that the calculated 2i//3/3 half-life is the same as its 
experimental value. The resulting gl^^^ is then almost the 
same one as with the old parametrization with the single gpp 
value. 

When the new parametrization of the particle-particle 
renormalization constants is used in the QRPA evaluation 
of the Oi^/3/3 nuclear matrix elements, a substantial reduction 
of the Fermi part, , is observed, while the Gamow- Teller 
and Tensor parts remain essentially unaffected. The full ma- 
trix elements M^^ are reduced by ~ 10 -20% as seen in Fig. 
[Sj We believe that such reduction, which also brings the ra- 
tio Xf closer to ~ —1/3, nearer to its value in the isospin 
conserving nuclear shell model values, is realistic, and should 
be used in the future application of the QRPA and its gen- 
eralizations. 
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TABLE II: Nuclear matrix elements for both decay modes with the old parametrization (grj^ ° = g^p ^ = Qpp) are compared to 
those with the new one, gjp"" 7^ ffJfT^- The adopted values of the parameter Qp^^ are also shown. The results for two values of the 
axial coupling constant qa arc displayed; the quenched value gA ~ 1-0 and the standard value gA = 1.27. The G-matrix elements of 
realistic Argonne V18 potential nucleon-nucleon potential are considered. The nuclear radius R = roA^^"^ with ro = 1.2fm is used. 



nucleus gA gpp " param. 2i^/3/3-decay NMEs Oi//3|8-decay NMEs xf 
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TABLE III: The same as in Table [TI] for realistic charge-dependent Bonn potential (CD-Bonn), 
nucleus qa gjp"" param. 2i'/3/3-decay NMEs Qvj3P-Aeca.y NMEs xf 
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1, 


.00 


0, 


.735 


old 


0.060 


0.145 


-1, 


.614 


3 


.035 


-0.373 


4.276 


2.655 


-0, 


.53 








0, 


.747 


new 


0.00 


0.145 


-1, 


.277 


2, 


.861 


-0.370 


3.769 


2.341 


-0, 


.45 




1, 


.27 


0, 


.753 


old 


0.055 


0.090 


-1, 


.583 


2, 


.640 


-0.351 


3.271 


3.271 


-0, 


.60 








0, 


.760 


new 


0.00 


0.090 


-1, 


.278 


2, 


.511 


-0.347 


2.957 


2.957 


-0, 


.51 




1, 


.00 


0, 


.770 


old 


0.097 


0.373 


-2, 


.883 


5, 


.532 


-0.615 


7.800 


4.8435 


-0, 


.52 








0, 


.769 


new 


0.00 


0.373 


-2, 


.366 


5, 


.526 


-0.606 


7.287 


4.525 


-0, 


.43 




1, 


.27 


0, 


.778 


old 


0.093 


0.232 


-2, 


.863 


4, 


.950 


-0.580 


6.148 


6.148 


-0, 


.58 








0, 


.776 


new 


0.00 


0.232 


-2, 


.367 


4, 


.950 


-0.571 


5.850 


5.850 


-0, 


.48 




1, 


.00 


0, 


.715 


old 


0.081 


0.423 


-2, 


.799 


6, 


.046 


-0.531 


8.314 


5.163 


-0, 


.46 








0, 


.713 


new 


0.00 


0.423 


-2, 


.288 


6, 


.052 


-0.521 


7.820 


4.856 


-0, 


.38 




1, 


.27 


0, 


.734 


old 


0.075 


0.263 


-2, 


.755 


5, 


.324 


-0.506 


6.529 


6.529 


-0, 


.52 








0, 


.732 


new 


0.00 


0.263 


-2, 


.290 


5, 


.330 


-0.497 


6.255 


6.255 


-0, 


.43 




1, 


.00 


0, 


.785 


old 


0.010 


0.206 


-1, 


.707 


3 


.942 


-0.253 


5.396 


3.351 


-0, 


.43 








0, 


.784 


new 


0.00 


0.206 


-1, 


.639 


3 


.940 


-0.251 


5.328 


3.308 


-0, 


.42 




1, 


.27 


0, 


.815 


old 


0.006 


0.128 


-1, 


.680 


3 


.564 


-0.241 


4.367 


4.367 


-0, 


.40 








0, 


.814 


new 


0.00 


0.128 


-1, 


.642 


3 


.563 


-0.240 


4.343 


4.343 


-0, 


.46 




1, 


.00 


0, 


.557 


old 


0.127 


0.200 


-1, 


.871 


4, 


.208 


-0.311 


5.768 


3.582 


-0, 


.44 








0, 


.555 


new 


0.00 


0.200 


-1, 


.057 


4, 


.206 


-0.305 


4.958 


3.079 


-0, 


.25 




1, 


.27 


0, 


.708 


old 


0.085 


0.00 


-1, 


.569 


2, 


.579 


-0.324 


3.230 


3.230 


-0, 


.61 








0, 


.707 


new 


0.00 


0.00 


-1, 


.062 


2, 


.575 


-0.321 


2.913 


2.913 


-0, 


.41 




1, 


.00 


0, 


.694 


old 


0.130 


0.0776 


-2, 


.673 


4, 


.902 


-0.580 


6.996 


4.344 


-0, 


.54 








0, 


.693 


new 


0.00 


0.0776 


-1, 


.850 


4, 


.887 


-0.572 


6.164 


3.828 


-0, 


.38 




1, 


.27 


0, 


.704 


old 


0.125 


0.0481 


-2, 


.641 


5, 


.582 


-0.549 


5.582 


5.582 


-0, 


.47 








0, 


.703 


new 


0.00 


0.0481 


-1, 


.851 


4, 


.476 


-0.542 


5.084 


5.084 


-0, 


.41 


130yg 


1, 


.00 


0, 


.698 


old 


0.102 


0.0545 


-2, 


.354 


4, 


.213 


-0.531 


6.036 


3.748 


-0, 


.56 








0, 


.697 


new 


0.00 


0.0545 


-1, 


.637 


4, 


.198 


-0.525 


5.310 


3.297 


-0, 


.39 




1, 


.27 


0, 


.707 


old 


0.098 


0.0339 


-2, 


.328 


3 


.867 


-0.502 


4.810 


4.810 


-0, 


.60 








0, 


.706 


new 


0.00 


0.0339 


-1, 


.637 


3 


.852 


-0.496 


4.373 


4.373 


-0, 


.42 




1, 


.00 


0, 


.664 


old 


0.110 


0.10 


-2, 


.268 


4, 


.532 


-0.486 


6.414 


3.983 


-0, 


.50 








0, 


.663 


new 


0.00 


0.10 


-1, 


.595 


4, 


.516 


-0.479 


5.632 


3.497 


-0, 


.35 




1, 


.27 


0, 


.712 


old 


0.092 


0.00 


-2, 


.231 


3 


.608 


-0.472 


4.522 


4.522 


-0, 


.62 








0, 


.712 


new 


0.00 


0.00 


-1, 


.599 


3 


.593 


-0.466 


4.119 


4.119 


-0, 


.44 




1, 


.00 


0, 


.657 


old 


0.0643 


0.0131 


-1, 


.300 


2, 


.395 


-0.269 


3.426 


2.127 


-0, 


.54 








0, 


.657 


new 


0.00 


0.0313 


-0, 


.858 


2, 


.385 


-0.268 


2.975 


1.847 


-0, 


.36 




1, 


.27 


0, 


.667 


old 


0.0619 


0.0194 


-1, 


.282 


2, 


.190 


-0.255 


2.735 


2.735 


-0, 


.59 








0, 


.667 


new 


0.00 


0.0194 


-0, 


.858 


2, 


.181 


-0.254 


2.460 


2.460 


-0, 


.39 
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TABLE IV: Ratio xf = Mp'/M^'^j. (see our definition of xf in tiie Eq. ^ ) in ISM [TD], QRPA-A, QRPA-B (present work, = 1.00 
and gA = 1-27 side by side), QRPA-A results are with the Argonne V18 potential, QRPA-B with the Bonn CD potential), IBM-2 [20], 
and QRPA-JyLa [21]. 

Nucleus X F 

ISM QRPA-A QRPA-B IBM QRPA-JyLa 





-0. 


,22 


-0 


.42, 


-0.51 


-0, 


.42, 


-0.51 


-0. 


.68 


-0.90" 


^«Ge 


-0. 


16 


-0 


.31, 


-0.34 


-0, 


.30, 


-0.34 


-0, 


,61 


-0.35 




-0. 


16 


-0, 


.33, 


-0.36 


-0 


.33, 


-0.36 


-0, 


.68 


-0.45 


'"Zr 






-0 


.45, 


-0.52 


-0, 


.45, 


-0.51 


-0, 


.10 


-0.69 


'""Mo 






-0 


.44, 


-0.49 


-0, 


.43, 


-0.48 


-0, 


,10 


-0.64 




-0. 


,25 


-0 


.39, 


-0.44 


-0, 


.38, 


-0.46 


-0, 


.05 


-0.61 




-0. 


30 


-0 


.43, 


-0.48 


-0, 


.42, 


-0.46 


-0, 


.10 


-0.45 




-0. 


,20 






-0.43 






-0.41 


-0, 


,56 


-0.68 




-0. 


,20 


-0 


.39, 


-0.43 


-0, 


.38, 


-0.41 


-0, 


,55 


-0.60 


130yg 


-0. 


,20 


-0 


.41, 


-0.44 


-0, 


.39, 


-0.42 


-0, 


,55 


-0.60 


"6Xe 


-0. 


,20 


-0 


.38, 


-0.41 


-0, 


.36, 


-0.39 


-0, 


,55 


-0.60 



" Ref.t22j 



